Magnetoelectric effect in multiferroic NdMn2O5
Bukhari Syed Hamad, Ahmad Javed
Department of Physics, Bahauddin Zakariya University, Multan 60800, Pakistan

 

† Corresponding author. E-mail: bukhari.hamad@gmail.com

Abstract

We have measured the dielectric constant for NdMn2O5 in an external magnetic field to map out the magnetoelectric phase diagram. The phase diagram corresponds well with the previously reported data of neutron diffraction and magnetic susceptibility. Our main finding is the observation of a dielectric anomaly in the low temperature phase with a strong magnetoelectric effect, which is attributed to the independent Nd ordering. Moreover, the absence of the dielectric anomaly in the paramagnetic phase is discussed, keeping in view the exchange interaction and its dependence on the rare-earth ionic radius.

1. Introduction

Magnetoelectric multiferroic compounds are a subject of significant consideration because they are characterized simultaneously by different magnetic and electric ferroic orders that are coupled. The magnetoelectric effect resulting from such strong coupling is stemming for exceptional materials, which opens a new horizon for the next generation applications in functional devices. In recent years, the magnetically frustrated multiferroic RMn2O5 family of oxides have attracted a great deal of attention due to the strong coupling between the magnetic and electric ferroic orders. [1, 2] The magnetoelectric effect revealed in RMn2O5 is due to a variety of phenomena, among them are magnetically induced polarization ( ) in HoMn2O5, [1] magnetically reduced in ErMn2O5, [1] magnetically flipped in TbMn2O5 [2] and GdMn2O5, [3] and magnetically flopped in TmMn2O5 [4] and YbMn2O5, [5] all arising in the low temperature phase. A reversible magnetocaloric effect has also been observed in HoMn2O5. [6] Moreover, the magnetoelectric coupling effect has also been studied in thin films. [79]

The main structure of all manganese RMn2O5 compounds crystallizes in the isostructural orthorhombic symmetry with space group Pbam. [10, 11] All members of RMn2O5 are characterized as the geometrically frustrated system because a slight deviation from their mean structure has been detected. [12] In a unit cell, there exist two manganese ions with different valances and positions: 4h for Mn (t e ) forming Mn O5 square pyramids and 4f for Mn (t e ) with corner sharing Mn O6 octahedra. These manganese ions together with oxygen coordination form a loop of five Mn ions, three Mn and two Mn . In each loop, there are three inequivalent magnetic superexchange interactions in the ab plane between the Mn and Mn ions, which are usually referred to as magnetic interactions J 3, J 4, and J 5. The J 3 and J 4 magnetic interactions correspond to the exchange interactions between the magnetic moments of Mn and Mn , and J 5 corresponds to the magnetic interaction between two Mn magnetic moments. Considering the ac plane, one can easily outline two Mn –Mn magnetic interactions, namely, J 1 and J 2. The J 1 is responsible for the magnetic interaction between adjacent Mn spins through oxygen at the layer and depends on the radius of the ion, which directly affects the z component (k z ) of the propagation vector ( ). [13, 14] Blake et al. have studied the crystal and magnetic structures of RMn2O5 ( , Ho, Dy) as a function of temperature by using neutron diffraction. They found that the spin arrangement of Mn ions is identical within the plane of the commensurate phase for each system. It was also suggested that a shift of the Mn cations gives rise to a net polarization along the axis due to the canted antiferroelectric structure. [15] Moreover, the radius of is essentially insensitive in the ab plane, while it determines the exchange interaction with adjacent Mn spins in the plane.

The key issue in the family of RMn2O5 is to elucidate the influence of on the magnetoelectric properties and its effect on the Mn magnetic ordering through J 1. It is important to note that RMn2O5 compounds with large size (La and Pr) are not ferroelectrics. [16] Interestingly, compounds with small size (Nd to Yb) all present a similar behavior, with the largest polarization in the commensurate (CM) phase, [2] except in the case of NdMn2O5. Very recently, it was observed that the multiferroic NdMn2O5 exhibiting the largest polarization along the b axis in the incommensurate (ICM) phase [17] differs from other RMn2O5. This makes NdMn2O5 an interesting compound to study because it lies between the nonferroelectric PrMn2O5 and ferroelectric SmMn2O5 and possesses a different magnetic origin for the ferroelectric polarization.

A spin–phonon interaction has been observed in the paramagnetic phase of RMn2O5 ( , Dy, and Bi) with the detection of Raman phonons [18] and infrared active phonons in DyMn2O5. [19] More interestingly, it was found that the phonon anomalies in the paramagnetic phase show a strong dependence on the ions. [18] Very recently, we observed an anomalous behavior in the dielectric constant slightly above the Néel temperature in GdMn2O5 [20] and DyMn2O5. [21] In this regard, we are motivated to conduct temperature dependent dielectric measurements in the paramagnetic phase of NdMn2O5 to reveal the existence of a dielectric response at high temperature in this compound. As Nd has a larger ionic size as compared to Gd and Dy , so we expect some interesting results that can shed light on the role of the ion in the family of RMn2O5.

Here, we report on the dielectric properties of polycrystalline NdMn2O5 in the temperature range T = 2–60 K with applied magnetic field H = 0–12 T. Based on our results, we provide a comprehensive explanation about the magnetoelectric transition in the NdMn2O5 system and also construct the magnetoelectric phase diagram. In addition, we will emphasize on the role of the ion in the paramagnetic phase.

2. Experimental setup

A polycrystalline sample of NdMn2O5 was prepared by using the sol–gel method. [20, 22] X-ray diffraction (XRD) measurement was performed with a Rigaku powder x-ray diffractometer using Cu radiation in the 2θ range from 20° to 80° by step scanning at 0.02° per step. The Rietveld refinement of the XRD pattern was performed by using the JANA2006 program. The pseudo-Voigt function was used to fit the diffraction peaks in the space group Pbam. The complex dielectric constant was measured in the temperature range T = 2–60 K with applied magnetic field H = 0–12 T at a frequency of 10 kHz by using a physical property measurement system (PPMS). For the dielectric measurements, the powder was pressed into a pellet of diameter 13 mm under the pressure of 30 kN by using the Paul-Otto Weber hydraulic press. The silver paste was used for forming a capacitor.

3. Results and discussion

The high quality single crystal of NdMn2O5 is difficult to grow by using the flux method. [23, 24] However, the polycrystalline NdMn2O5 has been synthesized by the chemical transport method under high oxygen pressure. [25] Here, we report the successful synthesis of polycrystalline NdMn2O5 by using the sol–gel method for the first time. Figure 1 presents the XRD pattern of polycrystalline NdMn2O5 at room temperature. The XRD pattern is analyzed by using the Rietveld refinement technique assuming an orthorhombic structure with Pbam space group (see inset of Fig. 1). A good agrement has been found between the refined and the observed XRD patterns as all XRD peaks are fitted, which indicates the single phase formation of the material and no trace of any impurity. The obtained unit cell parameters are Å, Å, and Å, which are consistent with the previously reported values. [12]

Fig. 1. (color online) Room temperature XRD pattern of polycrystalline NdMn2O5. Inset shows the Rietveld refinement of the XRD curve by using the JANA2006 program. The calculated crystallographic data are , , and .

The study of dielectric properties is an essential tool to identify the multiferroicity. Figure 2 shows the temperature variation of the dielectric constant during both heating and cooling with zero magnetic field. The dielectric constant increases sharply at ∼31 K, which is referred to as the Néel temperature. [17, 26] Just below , there is a sharp peak indicating a small thermal effect as observed previously. [17, 26] This sharp peak at K is attributed to the ferroelectric peak as observed in other RMn2O5. [27] It is worth mentioning that the thermal hysteresis of indicates the first order phase transition similar to the other ferroelectric members. A slight change in the dielectric curve observed at K may be associated with another weak dielectric anomaly.

Fig. 2. (color online) Dielectric constant ε of NdMn2O5 at zero magnetic field. The arrows indicate the transition temperatures between different magnetic phases.

If this anomalous behavior appeared as a sharp step, it might be attributed to excitation of electromagnons. [28, 29] Here, we suggest that the anomaly at K may not be related to any electromagnon excitation, rather it is related to a slight change in the magnetic propagation vector. [17, 26] The most interesting finding is the observation of the low temperature anomaly at ∼4 K, which is attributed to the independent Nd magnetic ordering as observed in other RMn2O5 ( Tb, [2] Gd, [20] Dy, [21] Er, [30] Tm, [31] and Yb [24]). Chattopadhyay et al. failed to observe the low temperature independent Nd ordering in their dielectric measurements. [17, 26] However, we have clearly observed Nd ordering in our dielectric measurements, which is consistent with the magnetic susceptibility and neutron diffraction measurements on NdMn2O5. [17, 26]

We have also measured the dielectric constant under a strong magnetic field in order to identify the magnetoelectric effect in NdMn2O5. Figure 3 demonstrates the variation in the dielectric constant under a magnetic field up to 12 T. Interestingly, the dielectric constant remains robust against the external magnetic field at the most prominent phase transitions ( , , and T X ), however the anomaly at smoothly vanishes with increasing magnetic field. The anomaly at shows an interesting magnetoelectric behavior and is quite similar to that observed for other members of the RMn2O5 family. [20, 21] Astonishingly, no magnetoelectric excitation is observed in the paramagnetic phase (above ) as observed in the case of GdMn2O5 [20] and DyMn2O5 [21] that gradually moves to higher temperature with increasing field. This reflects the role of the ions in the paramagnetic phase of RMn2O5, which is an unusual result and the origin of this discrepancy will be discussed later.

Fig. 3. (color online) Temperature dependent dielectric constant ε for NdMn2O5 at various magnetic fields up to 12 T.
3.1. Magnetoelectric phase diagram

In order to elucidate the magnetic field effect on NdMn2O5, we have constructed the magnetoelectric phase diagram, as shown in Fig. 4. The phase boundaries defined by , , and T X in the phase diagram remain robust under the external magnetic field. These phase boundaries correspond to the phase transitions PM/PE ICM1/PE at , ICM1/PE ICM1/FE at , and ICM1/FE ICM2/FE at T X , respectively (we have used the conventions for different phases (i.e., ICM1, ICM2, PE, FE) from Ref. [17]). At the lowest temperatures ( K), an additional structure is seen due to an ordered state of the Nd spins, which corresponds to another phase transition ICM2/FE CM+ICM2/FE.

Fig. 4. (color online) Magnetoelectric phase diagram of NdMn2O5.

The phase diagram for H = 0 exhibits similar phase sequences as observed through neutron diffraction. [17] Therefore, we suggest that fundamental physics behind the different phases is the same as that in the NdMn2O5 phase diagram as described in Ref. [17]. Above 3 T, the CM+ICM2/FE phase becomes unstable and merges with ICM2/FE as the field increases up to 12 T. A possible explanation for this unstable phase may be a reorientation of the Nd sublattice under the external magnetic field.

3.2. Discussion

The behavior of the dielectric constant is in line with the previous magnetic susceptibility and the neutron diffraction experiments on NdMn2O5. [17] The dielectric constant above 10 K is quite similar to the previously observed value and has shown no significant magnetoelectric effect. In most RMn2O5, the magnetoelectric effect appears in the low temperature phase due to different ions as observed here in the case of NdMn2O5. The magnetoelectric phase diagram clearly explains the magnetoelectric effect in the low temperature where the phase transition due to the Nd sublattice smoothly vanishes with increasing field and merges with the well-defined ferroelectric phase. Further investigation on the field dependent polarization should be carried out in the low temperature phase.

Flores et al. have revealed the existence of phonon anomalies in the paramagnetic phase just above for RMn2O5 ( , Eu, and Dy). They have found out that the phonon anomalies in the paramagnetic phase show strong dependencies with the ionic radii of the ions. [18] Thus it reflects that the paramagnetic phase in RMn2O5 is also sensitive to the particular ions as well as the low temperature phase. However, no detailed investigation has yet been made in the paramagnetic phase of RMn2O5. Recently, we have clearly observed a sizeable change in the dielectric anomaly in the paramagnetic phase (45–65 K) with strong field dependencies in GdMn2O5 [20] and DyMn2O5. [21] The absence of the dielectric anomaly in the paramagnetic phase (35–60 K) of NdMn2O5 is worth noticing and may be related to the bigger ion Nd (2.076 Å) as compared to Gd (1.056 Å). Two possible scenarios arise due to different ionic radii of . First, the bigger ionic radius of Nd shortening the bond length between adjacent Mn –Mn spins through oxygen at the Nd layer directly affects the k z component of the propagation vector as described in Ref. [14]. This increases the J 1 exchange interaction, which differs appreciably from the one observed in RMn2O5 with small ionic radius. Secondly, there is also a sizable effect in the ab plane that increases the J 3 while it decreases the J 4 exchange interaction due to the bigger ionic radius of Nd . [14] We assume that plays a significant role even in the paramagnetic phase due to a different ionic radius that changes the bond length and thus affects the superexchange interactions.

Currently, we are unable to speculate further about this discrepancy that arises in the paramagnetic phase. A detailed investigation of neutron diffraction may help to explain the role of the ions in the paramagnetic phase. Also, a study of spin-phonon coupling in NdMn2O5 through Raman and infrared active phonons with special emphasis in the paramagnetic phase would also be useful to confirm the role.

4. Conclusion

We have synthesized polycrystalline NdMn2O5 and performed dielectric measurements under a strong magnetic field. Our results unveil an anomaly in the low temperature phase, attributed to the Nd magnetic ordering. To summarize our results, a magnetoelectric phase diagram was constructed. We found that the main transitions defined by , , and T X in the phase diagram remain robust under the external magnetic field. However, a strong magnetoelectric effect is observed for the transition at that vanishes at fields higher than 3 T. Moreover, the absence of the dielectric anomaly in the paramagnetic phase may be due to the bigger ionic size of the Nd ion, which directly influences the J 1 exchange interaction and results in a different magnetic correlation in the paramagnetic phase as compared to the small ionic radius of in the family of RMn2O5.

Reference
[1] Higashiyama D Miyasaka S Tokura Y 2005 Phys. Rev. B 23 064421
[2] Hur N Park S Sharma P A Ahn J S Guha S Cheong S W 2004 Nature 429 392
[3] Lee N Vecchini C Choi Y J Chapon L C Bombardi A Radaelli P G Cheong S W 2013 Phys. Rev. Lett. 110 137203
[4] Fukunaga M Sakamoto Y Kimura H Noda Y Abe N Taniguchi K Arima T Wakimoto S Takeda M Kakurai K 2009 Phys. Rev. Lett. 103 077204
[5] Fukunaga M Sakamoto Y Kimura H Noda Y 2011 J. Phys. Soc. Jpn. 80 014705
[6] Ge H Zhang X Q Ke Y J Jin J L Liao Z X Cheng Z H 2013 Chin. Phys. B 22 057502
[7] Tang Z H Tang M H Lv X S Cai H Q Xiao Y G Cheng C P Zhou Y C He J 2013 J. App. Phys. 113 164106
[8] Tang Z Wang B Yang H Xu X Liu Y Sun D Xia L Zhan Q Chen B Tang M Zhou Y Wang J Li R 2014 App. Phys. Lett. 105 103504
[9] Tang Z Xiong Y Tang M Xiao Y Zhang X Yuan M Ouyang J Zhou Y 2014 J. Mater. Chem. C 2 1427
[10] Kagomiya I Kohn K Uchiyama T 2002 Ferroelectrics 280 131
[11] Alonso J Casais M Martínez-Lope M Martinez J Fernáandez-Díaz M 1997 J. Phys.: Condens. Matter 9 8515
[12] Balédent V Chattopadhyay S Fertey P Lepetit M Greenblatt M Wanklyn B Saouma F Jang J Foury-Leylekian P 2015 Phys. Rev. Lett. 114 117601
[13] Tachibana M Akiyama K Kawaji H Atake T 2005 Phys. Rev. B 72 224425
[14] Zobkalo I Gavrilov S Nyi N S Barilo S Shiryaev S 2014 J. Magn. Magn. Mater. 354 85
[15] Blake G Chapon L Radaelli P Park S Hur N Cheong S W Rodriguez-Carvajal J 2005 Phys. Rev. B 71 214402
[16] Doubrovsky C André G Gukasov A Auban-Senzier P Pasquier C Elkaim E Li M Greenblatt M Damay F Foury-Leylekian P 2012 Phys. Rev. B 86 174417
[17] Chattopadhyay S Balédent V Damay F Gukasov A Moshopoulou E Auban-Senzier P Pasquier C André G Porcher F Elkaim E 2016 Phys. Rev. B 93 104406
[18] García-Flores A F Granado E Martinho H Urbano R R Rettori C Golovenchits E I Sanina V A Oseroff S B Park S Cheong S W 2006 Phys. Rev. B 73 104411
[19] Cao J Vergara L Musfeldt J Litvinchuk A Wang Y Park S Cheong S W 2008 Phys. Rev. B 78 064307
[20] Bukhari S H Ahmad J 2016 Physica B 492 39
[21] Bukhari S H Ahmad J 2016 Physica B 503 179
[22] Zhang X W Lin T Xu J Xu L Chen K J 2012 Chin. Phys. B 21 018101
[23] Wanklyn B M 1972 J. Mater. Sci. 7 813
[24] Koyata Y Nakamura H Iwata N Inomata A Kohn K 1996 J. Phys. Soc. Jpn. 65 1383
[25] Euzen P Leone P Gueho C Palvadeau P 1993 Acta Crystallogr. C 49 1875
[26] Chattopadhyay S Balédent V Auban-Senzier P Pasquier C Doubrovsky C Greenblatt M Foury-Leylekian P 2015 Physica B 460 214
[27] Hur N Park S Sharma P A Guha S Cheong S W 2004 Phys. Rev. Lett. 93 107207
[28] Shuvaev A M Travkin V D Ivanov V Y Mukhin A A Pimenov A 2010 Phys. Rev. Lett. 104 097202
[29] Sushkov A B Aguilar R V Park S Cheong S W Drew H D 2007 Phys. Rev. Lett. 98 027202
[30] Kobayashi S Osawa T Kimura H Noda Y Kagomiya I Kohn K 2004 J. Phys. Soc. Jpn. 73 1031
[31] Kobayashi S Kimura H Noda Y Kohn K 2005 J. Phys. Soc. Jpn. 74 468